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We derive a closed set of equations for the kinetics and non-equilibrium dynamics of interacting Luttinger 
Liquids with cubic resonant interactions. In the presence of these interactions, the Luttinger phonons become 
dressed but still well defined quasi-particles, characterized by a life-time much larger then the inverse energy. 
This enables the separation of forward time dynamics and relative time dynamics into slow and fast dynamics 
and justifies the so-called Wigner approximation, which can be seen as a "local-time approximation” for the rela¬ 
tive dynamics. Applying field theoretical methods in the Keldysh framework, i.e. kinetic and Dyson-Schwinger 
equations, we derive a closed set of dynamic equations, describing the kinetics of normal and anomalous phonon 
densities, the phonon self-energy and vertex corrections for an arbitrary non-equilibrium initial state. In the 
limit of low phonon densities, the results from self-consistent Bom approximation are recaptured, including 
Andreev’s scaling solution for the quasi-particle life-time in a thermal state. As an application, we compute 
the relaxation of an excited state to its thermal equilibrium. While the intermediate time dynamics displays 
exponentially fast relaxation, the last stages of thermalization are governed by algebraic laws. This can be 
traced back to the importance of energy and momentum conservation at the longest times, which gives rise to 
dynamical slow modes. 
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I. INTRODUCTION 


The kinetics and non-equilibrium dynamics of low dimen¬ 
sional, interacting quantum systems is an outstand ing a nd fas¬ 
cinating challenge in quantum many-body physics On the 
one hand, it is strongly motivated by recent cold atom exper¬ 
iments performed on low entropy quantu m wires under out- 
of-equilibrium conditions! ^* l ^^ l "^^ l ^^ l ^^ l ^^ l ^^ l On the other hand, 
from a theoretical point of view, the study of non-equilibrium 
dynamics in integrable an d nonin tegrable systems is currently 
a field of growing interest. ^3131 This is triggered by the ques¬ 
tion, whether and - if answered affirmatively - in which 
specific w^ a one-dimensional quantum system is able to 
thermalize.cl 


An example of a one-dimensional integrable model is 
the linear Luttinger Liquid, which is the effective long- 
wavelength description of one-dimensional interact ing q uan- 
tum fluids, composed either of fermions or bosons.12313 Due 
to integrability, even if prepared in a non-equilibrium state, 
this model will never thermalize, since the number of excita¬ 
tions for each momentum mode ^ is a constant of motion^ 
However, as already pointed out by Andreev and Haldane,®^ 
there are non-zero corrections to the linear theory, which cer¬ 
tainly break integrability in the Luttinger model. They are ir¬ 
relevant in the sense of the renormalization group, and do not 
affect static observables. In contrast, they lead to a modifica¬ 
tion of dynamical, i.e. frequency resolved, correlation func¬ 
tions. These nonlinear corrections describe three-body scat¬ 
tering processes between phonons and are, due to their reso¬ 
nant nature, not straightforwardly approached theoretically. 

Apart from a wealth of numerical studies, 
based on matrix product state and Bethe ansatz cal¬ 
culations, several field the oretical approaches have been 
developedl ^^ l ^^ l ^ fi A seminal early study was car¬ 

ried out by Andreev, who used a self-consistent Born approx¬ 
imation to determine the phonon self-energies, establishing 


a universal phonon absorbtion rate ~ q’’ with exponent 
q - \ for a finite temperature system.®2122l a similar com¬ 
putation leads to an exponent 77 = 2 for the case of a zero 
temperature state,E3which has been verified by several numer¬ 
ical methods.l^ Recently, so-called nonlinear Luttinger Liq¬ 
uids have been introduced, which are designed to capture cor¬ 
rections to the linear Lutting er theory in the context of one 
dimensional fermions systems. EUSEUMl'pijese have been very 
sucessful in determining, for example, the power law diver¬ 
gences of the dynamic structure factor, or thermalization rates 
for near equilibrium systems.® Despite the lar ge number of 
analytic and semi-analytic works only few ap¬ 
proaches f or far fr om equilibrium dynamics have been devel¬ 
oped so far.®lll 2 ll 

The purpose of this article is to provide a quantitative 
description of the kinetics of an interacting Luttinger Liq¬ 
uid. Strong motivation comes from a recent surge of exper¬ 
iments with interacting one-dimensional bosons and fermions 
in ultracol d atom setu ps performed under out-of-equilibrium 
conditions 1|28|42|48 ]69] a one-dimensional quantum fluid 
is prepared in a true non-equilibrium state, and the experi¬ 
ment subsequently witnesses the time evolution of the sys¬ 
tem. At the present stage, the systems’ properties are still 
well described by linear Luttinger Liquid theory with, how¬ 
ever, a time dep endent, far from equilibrium distribution 
functionSo far, there are numerous theoretical works 
which describe these systems in terms of non-interacting Lut¬ 
tinger liquids alone, where there is only dephasing d ynamics 
and the evolution of the distribution function is absent.EISEIllI] 
This describes well the pre-thermalized regime in the shorter 
time dynamics. However, current experiments are steadily 
pushed to larger observation times, where the thermalization 
crossover caused by the residual (RG irrelevant) interactions 
in the Luttinger liquid should occur. What is therefore lack¬ 
ing at present in the theoretical literature, is a kinetic theory 
for the time evolution of the distribution function in these sys- 
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terns, which is able to track such a crossover and the associ¬ 
ated time scales. This is the aim of the present paper. 

On the technical side, achieving our goal amounts to obtain¬ 
ing theoretical control over the infrared divergences inherent 
to naive perturbation theory for the interacting Luttinger Liq¬ 
uid. While this issue was solved for the stationary equilibrium 
(or near equilibrium, in the sense of linear response) long time 
ago by Andreev and others,®3^we focus here on getting this 
problem under control in the kinetic equation, initialized with 
a general Gaussian non-equilibrium distribution function. 

In more detail, exploiting the resonant but subleading char¬ 
acter of the interactions in a one-dimensional interacting 
quantum fluid, we apply non-equilibrium diagrammatic the¬ 
ory to solve for the non-equilibrium dynamics of an inter¬ 
acting Luttinger Liquid. We show that due to the reso¬ 
nant but subleading nature of the interactions, vertex correc¬ 
tions are moderate for many physical realizations and conse¬ 
quently the non-equilibrium dynamics and self-energy can be 
solved within self-consistent Bom approximation. The self- 
consistency is however, crucial, and a perturbative Born ap¬ 
proach leads to infrared divergencies. The result is an effec¬ 
tively closed set of equations for the time-dependent phonon 
density and self-energy in the presence of resonant interac¬ 
tions. This approach, without considering the vertex correc¬ 
tion and restricted to equilibrium systems, has been discussed 
also for non-interacting one-dimensional dispersive fermionfl 
and in the context of the Coulomb drag effect for sufficiently 
low temperatures .1^ 

As a major result of the RG irrelevant but resonant in¬ 
teractions, the excitations remain well-defined but dressed 
phonons, with a life-time Tq much larger than their typical co¬ 
herent time-scale The dressed spectral function remains 
sharply peaked at the bare phonon energies Cq with width 
Tq^, such that the self-energy and distribution function of the 
phonons for frequencies sufficiently close to the on-shell fre¬ 
quency can be approximated by their on-shell value. This is 
referred to as the quasi-particle approximation. The long life¬ 
time of the dressed phonons results in a further simplification, 
as it implies that the forward time evolution of the system is 
much slower then the relative time evolution. This decoupling 
leads to a ’’local time approximation”, where an effectively 
stationary problem can be solved at each instant of time. We 
will quantify the validity of these approximations in terms of 
general but state dependent quantitative bounds below. 

Our estimate of the vertex correction further supports the 
validity of the self-consistent Born approximation for the time 
evolution. More precisely, for zero temperature states, the 
vertex correction vanishes identically, reproducing previous 
results.li^ Away from the ground state, the loop correction 
leads to a finite multiplicative renormalization of the vertex, 
which remains small for states close to thermal equilibrium. 
This implies that for typical translation invariant low entropy 
initial states, the equations governing the time evolution of the 
phonon occupation and the self-energy are effectively closed. 
It does not rule out, however, the possibility of significant ver¬ 
tex renormalization in general. 

The strength of this approach is the simplicity of the re¬ 
sulting final equations, which can directly be implemented 


and solved numerically. This provides a useful tool with a 
broad spectrum of applicability, ranging from tracking the 
thermalization process of non-integrable, weakly interacting 
Bose and Fermi ga ses to the study of interacting open system 
dynamics.It is suited for an initial state of a general 
Gaussian form with arbitrary, non-zero phonon densities, in¬ 
cluding diagonal and off-diagonal occupations. 

The structure of the kinetic equation and the equations for 
the self-energy and vertex correction reveal strong aspects of 
universal behavior. This concerns two key points: First, for 
an occupation function, which decreases sufficiently strongly 
in momentum space, they are independent of the short dis¬ 
tance cutoff of the Luttinger liquid, and thus of microscopic 
details inherited from even shorter length scales. This prop¬ 
erty emerges via a “bootstrap” mechanism within the self- 
consistent Bom approximation, as we elaborate on in the 
main text. Second, after a proper rescaling of timd^, all mi¬ 
croscopic parameters are completely eliminated from the dy¬ 
namic equations. Consequently and remarkably, the only mi¬ 
croscopic information that enters the dynamics is the initial 
phonon density at f = 0. 

At this point, we also mention three limitations of our ap¬ 
proach, which however are not directly physically relevant 
for the ultracold atomic systems we are targeting. First, as 
already briefly mentioned above, our setting is restricted to 
initial quantum states with Gaussian, although in general far- 
from-equilibrium correlations. Such c onditions hav e been 
discussed extensively in previous works * | 37 | 42|5Q|6 9] jj^gstly 
for the case of linear Luttinger Liquids generating the sub¬ 
sequent dynamics. Second, our method is not suited to de¬ 
scribe the asymptotic infrared behavior of a one-dimensional 
quantum fluid, where for the smallest moment^^ it has been 
shown that the elementary excitations are fermionJ^HEHSEMI 
with a strongly suppressed decay rate y ~ k^. Third, it 
doesn’t work for frequencies (or temperatures) above the Lut¬ 
tinger cutoff, where the dominant interaction is given by 
scattering between Luttinger phonons and mobile electronic 
impurities.l^Sl^However, the range of validity of our approach 
coincides perfectly with cold atom experiments, which con¬ 
sider momenta 10“^ A < p < A and temperatures k^T < Hu A, 
where u is the sound velocity and A the Luttinger cutoff. The 
Luttinger cutoff scale is given by the chemical potential for 
these experiments, A « sj2gnQ for weak interactions g, and hq 
the mean density. The infrared restriction of the momentum 
range results from typical trap sizes, with oscillator lengths 
roughly two orders of magnitude larger than 1/A. This dis¬ 
cards precision measurements of frequency resolved observ¬ 
ables in these extreme infrared asymptotic regimes, and thus 
is not of foremost interest for our study. 

The remainder of the paper is organized as follows. In 
Sec. [I^ we introduce and briefly discuss the action of the in¬ 
teracting Luttinger model in the phonon representation and 
in the Keldysh non-equilibrium framework. In Sec. Ill 


we derive the non-equilibrium fluctuation-dissipation rela¬ 
tion (FDR) and determine the phonon self-energies in self- 
consistent Born approximation for an arbitrary phonon distri¬ 
bution. Furthermore, we discuss the necessary approxima¬ 
tions and quantify their justification. Subsequently, taking 
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advantage of Sec. m we determine the kinetic equation for 
the phonon density in self-consistent Born approximation in 
Sec. IV In Sec.|V] we take into account non-zero off-diagonal 
(anomalous) phonon densities and show in which way the 
kinetic equation and the phonon self-energy are modified in 
their presence. Additionally, we derive a kinetic equation for 
the anomalous phonon densities. These results are applied 
in Sec.|V|]to determine the relaxation of an excited, thermal 
state back to thermal equilibrium. Furthermore, the analyt¬ 
ically obtained relaxation rate is compared to the numerical 
value, showing excellent agreement. Finally, in Sec. |VII[ we 
go beyond the self-consistent Born approximation and apply 
Dyson-Schwinger equations to take into account a non-zero 
vertex correction. We determine a closed set of equations for 
the kinetic equation, self-energy and vertex correction for ar¬ 
bitrary states, and discuss the effect of the latter. 


II. MODEL 


The action describing the interacting Luttinger model con¬ 
sists of two parts (we set ^ = 1) 

+ ( 2 . 1 ) 


Here, 5 tl is the w ell-known quadratic Tomonaga-Luttinger 
(TL) actioiPSmi 


“ 2;r 


- f 

Jx.t 


id_,(p) id,0) - uK {d,,9y - — {dxi/>y 
K 


( 2 . 2 ) 


where ^ = f_^ dt dx is the integral over space and time and 
(p — (p{x, t) and 0 - 9{x, t) are dimensionless, real fields. The 
non-linear part 5i„, is cubic in the fields and read^ 


^i X [^bc (d,9f (5,0) + {d,4>f\. (2.3) 


Starting from a microscopic derivation of the TL model as the 
effective long-wavelength description of interacting bosons 
or fermions in one dimensio n, the helds 9 ,0 represent local 
phase and density fluctuations.In this setting, we consider 
effective electron-electron (or boson-boson) interaction to be 
short ranged, i.e. of (5-function type.l^The non-linearity cor¬ 
responding to as well is of microscopic origin and is re¬ 
ferred to as band curvature. It originates from deviations from 
a perfectly linear dispersion of the microscopic particles. On 
the other hand, the term corresponding to is generated in an 
effective long-wavelength description, where the fast modes 
have been integrated out already, and describes effective three- 
particle interactions. 

The fields 0,0 are dimensionless, i.e. they have a canonical 
scaling dimension equal to zero. As a result, they do not scale 
when coarse graining to larger distances, i.e. when performing 
the rescaling 


X ^ lx, f —> ft, 


(2.4) 


where z = 1 is the dynamical exponent and I > 1. In contrast 


S Int ^ b 

2nl 


(2.5) 


under the rescaling ( |2.4[ ), such that the influence of 5i„, van¬ 
ishes on the longest wavelengths, i.e. it becomes irrele¬ 
vant in the renormalization group (RG) sense. Consequently, 
the static equilibrium properties of the interacting Luttinger 
model (Eq. ( |2.1[ )) are well described by the quadratic part of 
the action alone and the partition function can be approxi¬ 
mated by 


Z = 


^ lD[9,^]e'^ ^ J'£)[0,0]e'^TL (2.6) 


Here, J D[9,0 ] stands for the functional integral over the 
fields 9, 0. 

The Tomonaga-Luttinger action describes phonons with a 
dispersion Sg - u\q\ linear in the momentum \q\, propagat¬ 
ing with the speed of sound u. In the absence of 5i„,, these 
phonons are non-interacting and consequently the phonon 
density for a specific mode q is a conserved quantity. How¬ 
ever, although the phonon interaction is irrelevant in the RG 
sense, it contains resonant processes where two phonons prop¬ 
agating in the same direction and the same speed of sound can 
interact with each other for an infinite time span. This leads to 
a non-trivial modification of time-dependent, dynamical ob¬ 
servables compared to the case of non-interacting phonons. 
As pointed out in a seminal w ork b y Andree\0 (considering 
finite T) and more recent work,l^^l^ the presence of S i„, leads 
to a finite phonon lifetime 


Tq~\q\ '' with 


77 = I for r > 0 
,7 = 2 for r = 0 


(2.7) 


which, in equilibrium, is visible only in dynamical, i.e. fre¬ 
quency dependent quantities such as the dynamical structure 
factor 


S{q,u)^ f e'("''«")<(5,0),,,(5,0)o,o), (2.8) 

*Jx,t 


with 


{...)= ^ j£)[0,0]-e'". (2.9) 

For dynamical quantities it is therefore important to take into 
account the full action, Eq. instead of the reduced 

quadratic part only, as indicated in Eq. ( |2.6| l. 

Eor a true non-equilibrium situation, the phonon distribu¬ 
tion is not stationary, i.e. not a thermal or zero temperature 
distribution, but instead the phonon number tig for a given 
momentum mode q becomes time-dependent. The redistribu¬ 
tion of phonons between the different momentum modes with 
exact energy conservation is described by the resonant inter¬ 
actions in 5,„„ and in a non-equilibrium situation the action 
can not be reasonably reduced to the quadratic Tomonaga- 
Luttinger action, for which the phonon density is a constant 
of motion. 

Since we are interested in the non-equilibrium dynamics 
in the interacting TL model, we formulate the problem in 
a Keldysh path integral framework, which is able to treat 
both equilibrium and non-equilibrium dynamics on equal 
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footingP^I®^ We will now shortly introduce the canonical 
Bogoliubov transformation, which switches from the basis of 
real fields 6, (j) to the basis of complex fields a a. Those cor¬ 
respond to creation and annihilation operators in an operator 
picture!^ We close the model section by placing the action 
( |2.1| l on the Keldysh contour and briefly explaining the for¬ 
malism. 


A. Phonon basis 


In order to use a physically more appealing representation 
of the TL action, one commonly introduces a set of complex 
fields a a which represent the (bosonic) eigenmodes of the 
system, i.e. the discussed phonons. The corresponding Bo¬ 
goliubov transformation is 


0;c,r - ^0 + 2 ^ ^ ’ (2.10) 

1 

i f"' I '^ji, K. \ ^ 

^<Po - ^ J j sgn(^)e"'^^ (d^,, H- , (2.11) 

with abbreviations (fir , - (fi{x, t) and ^ ^ and the 

Fourier transformed phonon fields 




I 




( 2 . 12 ) 


The product ax,tax^t represents a phonon density and, there¬ 
fore, in the continuum limit, the fields Ox,,, Qxj are not dimen¬ 
sionless, in contrast to <fix,t,Qx,t, but scale as cixj ~ The 
quadratic part of the action transforms into 


S 


TL 


1 

27r 



{id, - u\q\)aq,,. 


(2.13) 


describing non-interacting phonons with a linear dispersion 
Cq = u\q\. The cubic part becomes 


5= r*- r Vq^p^p+q y/lqpip + q)\ 

2 ^ Jq,p,t 

f Cl-q-p^,Qq^tCipq 
^p+q,t^q,t^p,t + ^ + n 


with the vertex function 
qp 


^q,p,k — Kbc -J 2K 


kp qk 
\qp\ \kp\ \qk\ 


+ K, 


qp 


•C.), (2.14) 
(2.15) 


The interaction p.l4[ ) describes cubic phonon scattering pro¬ 
cesses with total momentum conservation. However, not all of 
the processes contained in S,„i are resonant, i.e. exactly energy 
conservi^ in the sense that ep+q — eq + ep. As explained by 
Andreev^ and pointed out above, the resonant processes lead 
to a divergence of the self-energy (and the kinetic equation, 
as we see later) in perturbation theory and are therefore the 
only relevant terms from a dynamical perspective. The term 
in Eq. ( |2.14[ ) describing the annihilation (creation) of three 
phonons can never be resonant. It will therefore play no role 


in our analysis and we will skip it from now on. For the resid¬ 
ual terms, resonance requires \p+q\ - |p|-i-|^|. For all momenta 
p, q fulfilling this condition, the vertex function takes on the 
value 


Vo = vi_ij 



(2.16) 


Consequently, instead of taking the full action 5,„i, it is suffi¬ 
cient to consider the reduced but resonant phonon interaction 


I y/\pq{p + q)\ [ap+q,taq,,apq + h.c.), (2.17) 

where the prime in ^ indicates that the integral runs only 
over momenta q, p which have the same sign. Together the 
quadratic action and the resonant phonon interaction describe 
the dynamics of the interacting Luttinger model in the phonon 
basis. 


+ (2.18) 


B. Keldysh action 


Non-equilibrium held theory is commonly performed in the 
Keldysh path integral framework, which is able to deal both 
with equilibrium and true non-equilibrium situations. To set 
up the Keldysh path integral, one first doubles the degrees of 
freedom in the theory by introducing plus and minus fields 
a+^q^,, a-^q^,, representin g for ward and backward time evolution 
on the Keldysh contour .I2E1 In this representation, the partition 
function is determined via 

Z = J (2.19) 

where S± is the phonon action p.l8| l with the replacements 
{ap,,dpq} {a±,p,t,d±^pj}. The+-representation contains re¬ 
dundancy, and a technically and physically more appealing 
representation is found by completing the transformation to 
the Keldysh representation, introducing classical and quantum 
fields according to 


flc = —= (a+ + a -), Gq - —= (a+ - a_). (2.20) 

V2 V2 


In the Keldysh representation, the quadratic action is 



with the bare inverse retarded/advanced propagators 

Dp,,,,' = S{t - t') (idp - u\p\ + 10+), (2.22) 

, = 6(t - f') {id,' - u\p\ - 10+) (2.23) 

and the Keldysh component of the inverse propagator 

D^,,, ^2iO^F{p,t,t'). (2.24) 
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Here, F{p, t, f) is the distribution function of the excitations 
and 0^ is the infinitesimal regularization for the quadratic 
theory.E] In an equilibrium, i.e. time-translational invariant 
situation, F(p, t, t') - F(p, t - f) and its Fourier transform is 
the bosonic distribution 

F{p, oj) = coth j = 2nB(w) + 1 (2.25) 


with the Bose function nsioj) - [e^ . The resonant 

interactions in the Keldysh representation take on the form 

^ f ^/|pk(k + p)\ [lal 4 al, 

VStt JpXt 

+ 44V) + ■ ^ 2 . 26 ) 

The bare response and correlation functions (retarded, ad¬ 
vanced and Keldysh Green’s functions) for the phonon de¬ 
grees of freedom are obtained according to 




- -i&(t - t') e 



, = /0(f' - t) 


^ -i{alXq,') - -{g^' ° 

0 g4 

'q,t,i' 


= -i (2nB{u\q\) -I- 1) (2.27) 


Here ... o ... stands for the convolution in the non-diagonal 
elements, i.e. the time index, but means multiplication in mo¬ 
mentum space. 

In the presence of interactions, the Green’s functions are 
modified by the emergence of non-zero self-energies 
which replaces the infinitesimal regularization. The corre¬ 
sponding formulas are 


G 

G 

G 


R 

q,r,r’ 

A 

K 


V /q,t,t' 

V 

- (g* o o G'^) 


(2.28) 


the bosonization procedure for interacting fermions or bosons 
out of equilibrium is in general not that simple, since initial 
conditions (even if quadratic in the microscopic fermionic or 
bosonic picture) will, in p rinciple, generate quantum vertices 
of arbitrary order,l2^E3llIl which have to be taken into account 
systematically for these cases. The aim of this work is to de¬ 
termine the self-energies and the distribution function 

Fiq, t, t') for a system that is driven out of equilibrium and 
evolves in time, for instance relaxing to an equilibrium state 
and approaching a bosonic distribution. 

This will be done in two parts. First, we show how one de¬ 
termines the self-energies from a given (non-) equilib¬ 

rium distribution function F{q, t, t'). To this end, we general¬ 
ize Andreev’s self-consistent Born approach to a non-thermal, 
non-equilibrium situation. Second, we use the kinetic equa¬ 
tion approach to determine the time-evolution of the distribu¬ 
tion function F in self-consistent Born approximation. Com¬ 
bining these two approaches allows us to determine the time- 
evolution of both the distribution function of the excitations 
and the self-energies, which play the role of finite lifetimes of 
the system’s excitations. 


III. SELF-ENERGIES 

The presence of the cubic, resonant interactions modi¬ 
fies the phonon response and correlation functions according 
to Eqs. ( |2.28[ ) by creating finite self-energies These 

self-energies are to leading order purely imaginary, leading to 
a finite decay rate of phonons or, in other words, to a finite 
phonon lifetime. We will now derive a method to determine 
these lifetimes for a non-equilibrium problem, where the dis¬ 
tribution function F{q, t, t') is time dependent and varies on 
time scales which are larger than the individual phonon life¬ 
times. To this end, we first derive the non-equilibrium version 
of a fluctuation-dissipation relation for the two-point response 
and correlation functions. 


A. Non-equilibrium fluctuation-dissipation relation 


where the infinitesimal factor (FF has been overwritten by the 
finite Keldysh self-energy Z^. The distribution function F in 
the presence of interactions is determined by the formula 

Gat a ={g‘^°F-FoG4 ■ (2.29) 

This setting corresponds to the physical situation, in which a 
system is initialized at time f = 0 in a Gaussian density ma¬ 
trix po (Gaussian in the bosonized language). It then evolves 
in time according to a Hamiltonian H, which constitutes the 
quadratic and cubic terms discussed in Eqs. ( |2.21| i and ( |2.26| l. 
In the Keldysh setting, this dynamics is expressed in terms 
of the retarded and advanced quadratic and cubic parts of 
the action. The initial density matrix enters the action in 
terms of pure quantum vertices. In the present case, i.e. for 
purely Gaussian initial conditions, the initial density matrix 
is completely captured in terms of the distribution function F 
and Eq. (|2.24[). However, as has been pointed out recently. 


Eluctuation-dissipation relations (EDR) relate the response 
(i.e. spectral) properties of the system encoded in G*^’^ to its 
correlations via the distribution function F. A particular ex¬ 
ample for such a relation is Eq. ( |2.29| l. Inverting this equation 
results in an EDR for the self-energies, 

Kut' = - -IJ^)oF-Fo{i4- z ^))^^ 

= -i{d, + dp)F,j,, + (z« O F - F O . (3.1) 

Eor a time-translational invariant system, the first term on the 
r.h.s. equals zero and due to the identity Z'^ = (z*) , Eq. ( |3.1| ) 
reduces to the well-known relation in frequency space 

E^, = -2/Im(z«4 F,,^. (3.2) 

This is consistent with our initial regularization of the 
quadratic sector for the case Z* = -iO^. 
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A useful representation for a two-time function F{q, t, t') is 
to choose Wigner coordinates in time, i.e. defi ning the for¬ 
ward time T = ^ and the relative time A, = f- f'PSThen one 
defines F{q,A,,T) = F(q,T + A,/2, t - Af/2) and its Fourier 
transform 


F(q, JdA,e‘^'‘^F(q,t,T). (3.3) 

Applying Wigner coordinates and the Wigner transformation 
p.3| l to the Keldysh self-energy in Eq. p.2| l leads to 

= -idrF,,^,r + o F - F C . (3.4) 

Eq. ^ is an exact expression for the Keldysh self-energy 
in Wigner representation. A complication arises since the 
Wigner transform of a convolution is not the product of the 
corresponding Wigner transforms, in contrast to the ordinary 
Eourier transform. In fact, one finds 




= s- 






drda^-da^dr 


q,co,T‘ 


(3.5) 


Without specific knowledge on the functional behavior of 2* 
and F, Eq. p.5| l is hard to evaluate explicitly. We will now 
briefly discuss a situation, with two particular approximations, 
which applies to the present model and for which the analytic 
evaluation of Eq. p.5|l is possible. 


1. Wigner approximation 


distribution at the cutoff temperature. We can thus apply the 
Wigner approximation to the EDR, resulting in 

= -idTFq^a,,T + 2/ Im (3.7) 

The validity of the EDR in Wigner approximation is a very 
important result. It is commonly used as the starting point 
for deriving a kinetic equation for the distribution function in 
arbitrary dimensions. However, in the present case, we will 
further simplify the EDR by making use of the fact that we 
are dealing with resonant interactions in one dimension. 


2 . Quasi-particle approximation and on-shell FDR 

The major effect of the non-linearities in the action is the 
emergence of finite phonon lifetimes due to resonant phonon- 
phonon scattering processes. The resonant character of this 
interaction - the fact that for two phonons travelling in the 
same direction momentum and energy conservation is ex¬ 
pressed by the identical 5-constraint - is the key property of 
one-dimensional systems with linear dispersion. The reso¬ 
nant contributions dominate the self-energy and the kinetic 
equation, while the non-resonant processes give only sublead¬ 
ing contributions to the lifetimes and the dispersion and have 
therefore already been eliminated on the basis of the action by 
using 5Res instead of 5i„, in Eq. p.l8[ ). 

The retarded self-energy is decomposed according to 

( 3 . 8 ) 


For the case of scale separation in the forward and rela¬ 
tive time, one can approximate the exponential in by the 
leading order terms. The product drdcj expresses the com¬ 
petition between relative time and forward time dynamics, 
it is smal l for slow forward time dynamics and fast relative 
dynamics .ES Comparing the zeroth order term with the first 
order term in an expansion of the exponential, we find the con¬ 
dition for approximating Eq. ( |3.5| l by the zeroth order term to 
be 


^oj^q,oj,T 

d F 

q,co,T 

q,(D,T 

J7 

q,ciJ,T 


If this condition is fulfilled, it allows for a separation into fast 
relative time dynamics and slow forward time dynamics. In 
general, it is a function of the Luttinger parameters u, K, the 
strength of the nonlinearity vq, as well as of the initial phonon 
distributions. Therefore, a general criterion for the applicabil¬ 
ity of the separation of time scales cannot be derived at this 
point. However, as has been shown in a different context,!^ 
for equilibrium or near to equilibrium initial states, it is ful¬ 
filled for temperatures k^T < uA smaller than the Luttinger 
cutoff. As it turns out, the right hand side of Eq. is a 
monotonic function of the quasi-particle distribution and a 
conservative criterion for the applicability of the Wigner ap¬ 
proximation is therefore riq < {e'^ - for all momenta q, 
i.e. the distribution should be smaller than the corresponding 


where cr^.^ is a positive, frequency independent function, 
which varies slowly in forward time r. Eor resonant inter¬ 
actions = 0 and consequently, the self-energy is fre¬ 

quency independent and purely imaginary. Non-resonant con¬ 
tributions in the interaction lead to 4 ^ 0, which however is 
generally strongly subleading compared to cr*. 

Eor the present model, the phonon interactions are RG- 
irrelevant and only their resonant character allows them to be 
of non-negligible influence. However, the effect of the inter¬ 
actions on the properties of the phonons will be small and 
subleading due to the RG-irrelevance. As a result, even in 
the presence of interactions, the phonons will have a lifetime 

Tp,, = -Im much larger than their associated coherent 

time-scale i.e. Consequently, the phonons 

remain well defined quasi-particles with a spectral function 
- i (G* - sharply peaked at the phonon energy 

(jj - Eor well defined quasi-particles, the self-energies 

YflAjK jjjg distribution function are evaluated on-shell, i.e. 
the frequencies are locked according to a; = u\q\, since fre¬ 
quencies CO + u\q\ do not contribute to the dynamics. Eor 
resonant interactions, on-shell evaluation is implied and 

sKitr = 0, (3.9) 

as stated above. This is consistent with the result of An¬ 
dreev’s and later works for equilibrium and holds true for the 
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/^K 


FIG. 1. Green’s functions in a diagrammatic representation. Full 
(dotted) lines represent classical (quantum) fields, while ingoing 
(outgoing) lines represent fields a (hermitian conjugate fields a). The 
time index has been omitted 


non-equilibrium case as well. The corresponding decomposi¬ 
tion for the Keldysh self-energy (with a convenient prefactor) 
reads 


= -2'0-^r + 


with ^ , = 0 


(3.10) 




. <' -rh.c. 


FIG. 2. Diagrammatic representation of see Eq. j2.26| l. In total, 
there exist six different vertices, three as depicted above and three 
corresponding hermitian conjugates. Each vertex has the prefactor 

Vip^q)- ^^^pqip + q)- 


and the on-shell self-energies 2*^ = = -2icr^ we 

can perform the frequency integration indicated in Fig. 1^ and 
find for momenta q > 0 


for resonant interactions. Inserting Eqs. p.9| l, p.l0| l in 
the non-equilibrium FDR p.7| l results in the on-shell, non¬ 
equilibrium FDR for resonant interactions 

CTg .^ = + crf,r (2n?,T + l) ■ (3.11) 

Here we have used the fact that for well defined quasi¬ 
particles, the on-shell distribution function F^,tj=u| 9 i,r = 2n^ ^ -H 
1 equals the time-dependent phonon density at momentum 

q- 

Eq. p. 1 l| i is the final form of the non-equilibrium EDR that 
we will use to set up the kinetic equation in the following sec¬ 
tion and to determine the Keldysh self-energy cr^ for a system, 
for which the time-dependent phonon density is known. 
Eor the latter case, the only unknown quantity is the retarded 
self-energy cr* and we can now set up the diagrammatic com¬ 
putation of cr® in the Keldysh non-equilibrium framework. 


B. Quasi-particle lifetimes in self-consistent Born 
approximation 

In this section, we perform the self-consistent Born approx¬ 
imation in a diagrammatic representation to obtain the self¬ 
energies cr*.^. This amounts to an infinite resumation over a 
certain class of diagrams and cures the divergence of the self¬ 
energy occuri'ing in the perturbative diagrammatic approach. 
In Sec. |VII[ we demonstrate that the self-energy can be de¬ 
termined exactly in a one-loop computation using Dyson- 
Schwinger equations and show that the deviation from the 
self-consistent Born approximation is negligible for many ini¬ 
tial states. 

In a diagrammatic approach, classical (quantum) fields are 
represented by a full (dotted) line, while ingoing lines rep¬ 
resent fields a and outgoing lines their complex conjugate a. 
This leads to a representation of Green’s functions in terms 
of diagrams as indicated in Eig. and vertices as depicted 
in Eig.|^ The retarded self-energy in self-consistent Born ap¬ 
proximation can directly be derived by common diagrammatic 
rules and is depicted in Eig.[^ Inserting the Keldysh Green’s 
function 


« _ 2 f f 

r pq(p + q)o-p ] 

Jo<p cr«(cr« 


/ 


pq{p - q)o-f, 


q<P cr« cr« + o-‘ 


.R ) 

p-qj 


(3.13) 


Since the self-energies must be invariant under the transfor¬ 
mation p —> -p, one can further simplify Eq. p.l3| l, ending 
up with 




=^0 r 

JOKp 


qp(q-p) , qpip + q) 


(cr« H-o-f 


q-p 


O-p+O-p+q/ 


(3.14) 


Einally, cr^ can be replaced using the EDR p.l l| l, which leads 
to 



+ 2.Tlp + 1 


qpiq - p) 


o-p + o-q-p 


qpip -b q) '' 

(3.15) 


For a given, time dependent distribution function 
is the only unknown function in this equation and has to be 
determined self-consistently. For a general time-dependent 
function ng^j this has to be done by iterating Eq. p.l5| l nu¬ 
merically until a self-consistent solution has been found. Eor 
the particular case for which ng^r shows scaling behavior in 
a sufficiently large momentum window, one can determine a 
scaling solution for the self-e nerg y as well and extract the cor¬ 
responding scaling exponent .*^1^ We will now briefly discuss 
the latter case and determine possible scaling solutions for the 
self-energy, and close the section with a discussion on univer¬ 
sal aspects of the scaling solution. 


1. Scaling solution for the self-energy 

For the case when the density Ug^ is a scaling function, i.e. 

(3.16) 

it is easy to show that also will be a scaling function 


^q,(jj ~ ^q,(jj^q,(jj^q,(jj 


( 3 . 12 ) 


= ( 3 - 17 ) 
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lP 9 (<l-p)i Gp Gq-p 

^ N(g + p)| G^+p 

^ \pq{q + p)\ GP+q 


FIG. 3. The retarded self-energy is the sum of three distinct diagrams 
as shown above. The notation is 2 = [a),q], P = {v,p] with p,q 
being momenta and a>, v the corresponding frequencies. The mo¬ 
mentum integral runs only over momenta that fulfill the resonance 
condition, which requires a) = u\q\. 


where the exponent rj^ and prefactor yf solely depend on the 
scaling behavior of i.e. on ?]„ and a^. In order to show 
this, we introduce the rescaled self-energy it* = ct^/vq and 
time T = t/vq, leading to 



FIG. 4. Integrand in the self-energy integral = £ ii,,i^(x) 

as a function of x = pjq, where q is the external momentum. The 
integrand j{x) elucidates the contribution to a£ from different mo¬ 
menta and shows clearly that the self-energy at momentum q is dom¬ 
inated by the behavior of cr^ for p < q. For physically relevant distri¬ 
bution functions, with the property —> 0 for ^ > A, the exponent 

—> 2 for X —> oo, which leads to a decay of the integrand ~ x“^ 
and a well-defined self-energy integral. 


2. Finite temperature states: For » 1 and 

q > maxfj^^, only the factor proportional to 

contributes and we have 

77« = 2+^, ^ (3.21) 



+ 2np 



qp(q - p) 


(T^ + (T 
^ p ^ ^ q-p 


+ 


qp{p + q) ' 

(3.18) 


Next, we insert Eqs. ( |3.16| ), ( |3.17| i into ( |3.18| l and perform the 
transformation p —> qx, yielding 


5 f /= 4 /-^ r 

yx Jo<x 


x(l - x) 


x(l -H x) 


X''" -H 11 - xl''" X''" H- |1 H-xl't", 

(3.19) 


X 11 H- lofiqx)’’" + ^^^(qx)’’" 

n 


The exponent 77 * is bounded from below and from above ac¬ 
cording to 1 < 77 * < 2. Here, 1 < 77 * results from the fact that 
the interaction is RG irrelevant and the self-energy correction 
can only be subleading compared to the dispersion - u\q\. 
The case 77 * = 2, i.e. diffusive scaling of the lifetimes is only 
reached in the zero temperature situation, otherwise one ex¬ 
pects superdiffusive behavior due to finite phonon densities 
with 77 * <2. As a result, Eq. p.l9| l distinguishes between 
three regimes: 

1. Low temperature states (T a; 0): Eor Uq ^ \ and 

drtiq ^ 0 the only term in the second row of Eq. p.l9[ ) with 
a relevant contribution is the constant unity and 


77« = 2, ^= 7744 =^, (3.20) 


Eor a finite temperature state, 77 ^ = ^ as \q\ 0, such that 



0.789 



(3.22) 


3. Non-equilibrium states: Eor 7 ?^ » 1 and 
q < mini , .. /), the dominant contribution stems from 

^ 'ZOrr/i \j2a,yR>' _ 

the time derivative in Eq. P.191 l and one finds 

77 * - ^ (3.23) 


The integral factor is defined as 

r / x'^'41 - x ) x'^'41 + zc)') 
Jo<xU'' + |1 - x|'i ’’’ X'l H- |1 H- X|'> j ■ 


(3.24) 


This gives rise to the possibility of generic non-equilibrium 
scaling behavior for out-of-equilibrium states, with an expo¬ 
nent 77 *, which is different from the zero and finite tempera¬ 
ture cases. A possible scenario of such scaling is discussed 
in Ref.for a driven system of interacting bosons with a non¬ 
equilibrium distribution with 77 „ = 1 in the infrared. 

The results discussed above include the case of finite and 
zero temperature equilibrium. Eor the latter riq - 0 and there¬ 
fore 77 * = 2 ,yf = Vo yjj, while for finite temperature ^ 

and ^^nq - 0 and consequently 77 * = 3/2, y^ - 0.789vo 
These are the mentioned results o btained for the zero and fi¬ 
nite temperature equilibrium cases. 


where the factor //j ,, will be defined in the following. 
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+ 




— ^ 9)1 i^P^Q-P + ^P^Q-p) 


+ 2 




]pip^'‘ 


FIG. 5. Diagrammatic representation of the Keldysh self-energy cr^ 
in self-consistent Born approximation. The index Q = ± ia^, q] 

represents the on-shell frequency and momentum relation (- 1 - for re¬ 
tarded, - for advanced Green’s functions), while P = lv,p] is the 
inner loop frequency and momentum and . Integration over 

resonant processes only is implied. 


2. Insensitivity of the self-energy to UV-behavior 


The insensitivity of the above results to the behavior of the 
model in the ultraviolet (UV) regime and therefore the non- 
universal properties of the system is guaranteed by the struc¬ 
ture of the vertex and the diagrams in Fig.|^ The self-energies 
cr^ are dominated by loop momenta p < q below the external 
momentum q. Therefore the non-universal behavior in the UV 
does not enter the self-energies. This is emphasized in Fig.j^ 
where the integrand 


|1 - x\‘' H- (1 -I- -H X‘' 


(3.25) 


in /,| ,2 = X - plq,is plotted for the case of thermal 

and zero temperature scaling. It is evident that for ingoing mo¬ 
menta q only momenta p < q contribute, which show the same 
scaling behavior in the distribution function. Therefore the 
scaling solutions for the self-energy are robust against modifi¬ 
cations of the density and the model itself when approaching 
the UV. 


IV. KINETIC EQUATION FOR THE PHONON DENSITY 

In this section, we will derive a kinetic equation for the 
phonon density in the case of resonant interactions. A kinetic 
equation describes the time-evolution of the distribution func¬ 
tion gen erated by the Keldysh and retarded/advanced self- 
energies,l 232 ll 2 l which can then often be evaluated in the per¬ 
turbative Born approximation. In the present case, due to the 
fact that the interactions are resonant, the kinetic equation per¬ 
turbation theory diverges and we have to use a self-consistent 
Born approach. Therefore, we evaluate the self-energy dia¬ 
grams with full Green’s functions, as in the previous section, 
which leads to an effective vertex correction for the kinetic 
equation and a time-evolution linear in the interaction param¬ 
eter Vo. 


The time-evolution of Uq^ is determined by the solution of 
the on-shell FDR in Eq. ( |3.1 l| i, which after rearrangement 
reads 


djHq = a-q .^ - (iriq^r + l) (4.1) 

Again, the retarded self-energy is determined via a diagram¬ 
matic approach, where the corresponding diagrams are shown 
in Fig. 1^ However, in contrast to the previous section, we 
also use a diagrammatic approach to determine the Keldysh 
self-energy cr^. As a consequence, we derive a non-linear dif¬ 
ferential equation for the distribution function. The diagram¬ 
matic representation of the Keldysh self-energy is depicted in 

Fig|5] 

For the full Green’s functions, we use the results from the 
previous section 


= 27r (w - u\q\ ± ‘ (4.2) 


and the relation 




= G 


* F 

q,OJ,T^ q,lP,T 




(2«?,t + l) 

(w - u\q\f -H [o-lrf 

(4.3) 


In Eq. ( |4.3| l, the first equality holds in Wigner approximation, 
while the second equality results from the quasi-particle ap¬ 
proximation, both discussed in the previous section. 

The frequency integration in the diagrammatic representa¬ 
tion can be performed analytically and yields the kinetic equa¬ 
tion (omitting time index) 


5i-n^-2vQ I {tipTiq-p nq[lrip-y riq-p^ 

I'o f («9 + Mp + l) - nqtip) (4.4) 


+4v" 


After transforming to dimensionless variables cr^ = vocr*, t = 
—, we finally arrive at 

Vn ’ 


dftlq- j yfSX {uptlq-p llq^l -y lip -y tlq^pjj 

Jo<p<q V V // 

+ r J^{npXn 


l)-nqnp). (4.5) 


This is the kinetic equation for the phonon density in self- 
consistent Born approximation. 

Comparing this equation to Eq. p.l8| l, one finds that the 
rescaled self-energy and therefore also the kinetic equa¬ 
tion for a rescaled time r = vqt only depend on the phonon 
distribution riq^f and is independent of all possible microscopic 
details that may enter u, K, vq in the model. As a result, the 
dynamics in the rescaled variables is very generic for an inter¬ 
acting Luttinger Liquid and only depends on the initial distri¬ 
bution function with which the system is initialized. 

The typical time-scale for the kinetic equation in the orig¬ 
inal variables is T,yp = i.e. linear in the vertex vq. This 
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is a non-perturbative effect resulting from the resonant inter¬ 
actions. Since two vertices enter the one-loop diagrams, one 
might naively (or in perturbation theory) expect that the typi¬ 
cal time scale is given by the square of the non-linearity. How¬ 
ever, this is normalized by the self-energies, which are propor¬ 
tional to Vo and required to regularize the kinetic equation. 

In the kinetic equation in ( |4.5[ ), still the self-energies o-* oc¬ 
cur. In principle, one could again replace them by a diagram¬ 
matic expression, which would give rise to an infinite hierar¬ 
chy. However instead of doing so, we use an iterative process 
in which, for a given time t, we determine self-consistently 
the self-energy cr^^. This result is then used to determine the 
r.h.s. of the kinetic equation, to subsequently compute the 
distribution function in the next time step r -H dr. This pro¬ 
cedure is illustrated in Fig. and allows us to compute the 
non-equilibrium dynamics of an interacting Luttinger Liquid, 
which may be initialized in a non-thermal state. 

We will now close this section by discussing two limiting 
cases where analytic results become available. First, the ki¬ 
netic equation for small external momenta q, and second the 
kinetic equation for a phonon distribution Uq j - riqj + 5nqj 
close to an equilibrium distribution Uq j. 


A. Kinetic equation for small external momenta 


For small external momenta, the kinetic equation ( |4.5| l can 
be brought into an even simpler form, explicitly revealing the 
evolution of for small q. In this case, the first integral in 
( |4.5| l covers only a very small momentum region and is pro¬ 
portional to q'^^^ , 1 < 77 * < 2. As a result, it is negligible 
compared to the second integral. On the other hand, in the 
second integral summations including q can be replaced ac¬ 


cording io p + q : 
then simplifies to 




+ 


CTp. The kinetic equation 


df^q,f — 1^1 f* "t" ^/j.f} — (4.6) 

Jo<p cr« - ^ ' 

where is a time dependent but momentum independent 
functional. The phonon density becomes 


q<s:l 


= nqj=o + \q\ 


f 


I, 


(4.7) 


0<f<f 


for sufficiently small momenta. Crucially, the change is linear 
in momentum q and for q - Q the density is time independent, 
and for all times pinned to its initial value. This is an exact re¬ 
sult, which can be traced back to particle number conservation 
in the underlying microscopic model.l^ 


B. Relaxation close to equilibrium 



next 
time step'' 


' d-j-'riq^T 2^ 


' kinetic 
equation 


FIG. 6. Schematic illustration of the iteration process to determine 
the time-dependent phonon density iiq^r- For a given time t, the self¬ 
energy is determined via the self-consistent Bom approximation 
according to Eq. ( |3.18k Subsequently, the time-derivative of riq.^ 
is computed via the kinetic equation ( |4.5k Using a Runge-Kutta 
solver for numerical differential equations, the density n^ r+sr can be 
computed and used as the starting point for the next iteration. 


for arbitrary temperature T. Sometimes one is interested in 
the relaxation of the distribution function close to equilibrium 


nqj = n^{.u\q\, T) -t- Snq^, (4.9) 


where the variation 6n Hb- For small momenta u\q\ T, 
nj^{u\q\, T) = ^ and we can expand the kinetic equation in the 
variation 6n. The zeroth order part solves the kinetic equation, 
and the leading order contribution stems from the first order 
of the expansion. After eliminating negligible terms, it reads 


d-^Suq = - 


ITq^dtiq r 

Jo- 

/ 


\J0<p<q <5-* + 0-« H- (T* 


q-p 

\ 


o<pd-^ + d-^ + o^+p/ 


(4.10) 


For a distribution close to thermal equilibrium, the self¬ 
energy will take on the thermal form and 


dfSuq X, —a-T 6nq 



— Q^T;O^Snq^ 


(4.11) 


where cf* is the thermal self-energy and ctt ~ 1.1 is a univer¬ 
sal number. This result can also be seen as an expansion of 
Eq. (|4.1|i in Stiq, which is 


dfdriq = -26nqd^ + 


dSn, 


l^d-q - (2ns + 1) 




6nq. 

(4.12) 


The second term thus gives a correction to the simple factor 
of 2 in Eq. ( |4.1 l| l. We thus obtain a nonperturbative estimate 
for the relaxation time of the interacting Luttinger Liquid 


A stationary solution of the kinetic equation = 0 is 

given by the Bose distribution function 


0.868 I u _3 


(4.13) 


nq,r ^ ns(u\q\,T) = - 1 ) 


reflecting the very slow asymptotic approach to equilibrium of 
the long-wavelength modes. It is very similar to the lifetime 
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of a single thermal phonon (cf. Eq. p.22| i), only modified 
by the prefactor a-^ - 1.1. This modification arises due to 
in-scattering processes of excitations p q scattering into 
the mode q. Since the main relaxation process is caused by 
out-scattering processes, the above result is quite intuitive and 
supports the statement that relative time dynamics ~ are 
fast compared to forward time dynamics ~ t^. 


V. KINETIC EQUATION AND DIAGRAMS IN PRESENCE 
OE ANOMALOUS DENSITIES 


In this section, we consider the effect of non-zero anoma¬ 
lous (off-diagonal) phonon density, i.e. agci-g 0 on the di¬ 
agonal kinetic equation and self-energy, and derive the kinetic 
equation for the anomalous density. Off-diagonal phonon den¬ 
sities can be populated due to external perturbations, as for 
instance density m odula tion due to a Bragg beam or a global 
interaction quench,f3Sl and their impact on the kinetics and 
non-equilibrium dynamics is therefore non-negligible. Sum¬ 
marizing the results of this section, due to the structure of 
the resonant interactions, the kinetic equation for the diago¬ 
nal phonon density and the diagonal retarded/advanced self¬ 
energies are not modified in the presence of anomalous densi¬ 
ties and remain unchanged compared to Eqs. ( |0] l and ( [TT5 ] i. 
In contrast, the kinetic equation for the anomalous densities is 
fed by the normal occupations, cf. Eq. ( |5.25| ). 

Eor a generic equilibrium situation and for certain real¬ 
izations of systems brought out of equilibrium, the anoma¬ 
lous (off-diagonal) phonon density and consequently the cor¬ 
responding response and correlation functions are zero, i.e. 

= (5.1) 

where a, a' - c,q represent classical or quantum indices. 
However, when a system is driven out of equilibrium, it is 
possible to populate off-diagonal terms. A simple situation 
for which this is indeed the case is an interaction quench in a 
one dimensional quantum fluid, where the off-diagonal densi¬ 
ties are non-zero after the quench and lead to non-equilibrium 
dynamics even in the absence of phonon-phonon scattering 
processes 

In order to deal with the situation of anomalous densities, 
we first have to modify the EDR accordingly. The Keldysh 
Green’s function in the presence of off-diagonal terms, ex¬ 
pressed in Nambu space, is 


G 


K 


( 


V 


Sq,t,t' 



q,tj' 

, 

o-q,t',t 


( KKr> \ 

(5.2) 


Eor the quadratic theory in the absence of phonon scattering, 
i.e. 5 = 5 tl only, the Keldysh Green’s function can be evalu¬ 
ated explicitly. In a operator representation, it reads 


G 


K 

q,l,t' 



<{V’' 


(5.3) 


with the anti-commutator {■, ■). In Wigner representation it is 

K _ Sico- €q)(2nq +1) d(w)2m^e~'^^’’' \ 

' \ (5(w)2m;e'2"-'^ 6(cj + eq)(2n^q -H 1) j ’ 

(5.4) 


where nig is the anomalous phonon density (nig — \{aga-g}\ in 
an operator representation in terms of annihilation operators 
a). In the quadratic theory, both the normal and the anomalous 
densities are constants of motion. 

The Keldysh Green’s function in ( |5.4| ) has two essential 
drawbacks. Eor the case of non-zero but slowly (compared to 
eg) varying anomalous density irig, the off-diagonal terms of 
G^ are not slow but oscillate with the fastest scale in the prob¬ 
lem, i.e. f,yp = 27 ' The Wigner approximation is therefore 
not applicable for the off-diagonal terms of G^. Eurthermore, 
the Keldysh Green’s function in frequency representation is 
peaked at three different frequencies, a> - (e^,0,-e^). Both 
prevents an EDR in the form of Eq. p.7[ ) to exist in this repre¬ 
sentation. 

In order to cure this problem, we switch to a rotating frame 
by introducing the fields Ogj = ag ,e“^‘i', , = ag ,e~“^‘i\ which 

modifies the quadratic action according to 



0 id, -t- iO^ 
id, - iO^ 2i0^ coth 


(5.5) 


The resonant phonon interaction is invariant under the trans¬ 
formation 

[^^k+pXk.t(^%r 

+^Lp,t + h-c.]. (5.6) 

since the phase = 1 in the case of resonance, i.e. 

for \k + p\ - |k| -I- \p\. The corresponding correlation function 
in Nambu space and Wigner coordinates after the rotation is 


Gg,0J,T = 


2na -1- 1 2nia 


2m* 


2n^g + 1 


(5.7) 


while the bare retarded Green’s function reads 
1 o \ 1 


G* — I w+(0+ 

I Q 


-i0+ 


= cr. 


0) + /0+ ’ 


(5.8) 


where cr^ is the Pauli matrix. Respecting the symplectic struc¬ 
ture in bosonic Nambu space, the EDR in the presence of off- 
diagonal densities is 




oo-^oF-Fo(t^o 


G^) = -iS(io)Fg 

'q,cj,T ^ 


(5.9) 


Here F is the physical distribution function for the phonons, 
with the on-shell value 

P —I '^^q,T 1 \ 

2m;, 2n_,., + lj' 


(5.10) 
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9q 


h^-Q 
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atg ag 


S{zpP + zp'P') 

an , / \Q^n' 7 , 

right vertex; 

S{zQ^q' + ZK'k' + zp'p') 
i{zQ,\q'\ + ZK'\k'\ + zp,\p'\) 


left vertex: c q 

HzQq + ZKk + zpp) “ff.zK ^K',z^, 

S(zQ\q\ +ZK\k\ +zp\p\) S(zkK + ZK'f{:') 


FIG. 7. Diagrammatic representation of the diagonal and off- 
diagonal retarded/advanced and Keldysh Green’s functions. In the 
presence of anomalous densities, the off-diagonal Keldysh compo¬ 
nents will be generally non-zero, t- 0. The retarded and advanced 
Green’s functions only differ from zero if a finite off-diagonal self¬ 
energy is present, cf. Eq. l|5.12[(. 


Both the transformation to a rotating frame by switching from 
{aq,aq] to {aq,aq] as well as the symplectic factors cr^ in 
Eq. ( |5.9| ) are necessary modifications in order to obtain an 
FDR with a physically relevant distribution function F. The 
latter should consist of diagonal and anomalous densities that 
are slowly varying in time and reproduce the matrix structure 
ofG^. 

Inversion of Eq. (|5.9|l yields the kinetic equation 


^djFq (p T — kT 


R -T' 


J Fq,til,T^q,ip,T^Z 


(5.11) 


In the absence of off-diagonal terms in both self-energies 
and the distribution function, this equation reduces to the or¬ 
dinary kinetic equation discussed in the previous section. We 
will now set up the diagrammatic computation of the self¬ 
energies in order to derive the kinetic equation in the 

presence of anomalous phonon densities. 


A. Diagrammatics for off-diagonal terms 


In the presence of off-diagonal densities, one can no longer 
generally exclude non-zero off-diagonal self-energies and 
consequently non-zero off-diagonal Green’s functions from 
the action. In this section we set up the diagrammatic com¬ 
putation of the self-energies in Nambu space. We obtain two 
key results, which crucially rely on the resonant character of 
the interaction. First, in the retarded/advanced sector, the off- 
diagonal self-energies are exactly zero even in the presence 
of anomalous densities. Second, we compute the off-diagonal 
self-energies in the Keldysh sector, which are different from 
zero. 

The retarded Green’s function in Nambu space is (we use a 
general index Q - {q,<jj,T), -Q - {-q,-a),T) for this para¬ 
graph) 


■R -I 4 




-Q ^-Q 








-P 

Q 

-OJ - 


{CO - E«)(m + E^g) + 


-p 


+ CO 

Q 


(5.12) 



FIG. 8. Diagrammatic representation of all possible one-loop di¬ 
agrams contributing to the retarded sector of the self-energy. The 
additional index Zq,p,k = ±1 represents ingoing lines {aQ,p,K fields) 
for z = 1 and outgoing lines (Hq^pk fields) for z = -1. The d- 
constraints stem from the momentum/frequency conservation in the 
Green’s function and the momentum/frequency conservation in the 
two vertices. Two additional constraints, one for each vertex, are 
caused by the resonance condition. This results in six crucial con¬ 
straints discussed in the main text. 


/Zg = -iXo-^gOTg). The diagrammatic representation for the 
retarded and Keldysh Green’s functions in Nambu space is 
depicted in Fig. |7] 

With this at hand, we can set up diagrammatic rules in 
Nambu space, as we do in the following. To this end, we start 
from the most general diagram contributing to the retarded 
sector of the Green’s function, which is shown in Fig. Here 
we replaced the fields ag ^ trgj and dg —» ag-i in or¬ 
der to find a generalized representation of diagrams in Nambu 
space. Exploiting frequency and momentum conservation in 
the Green’s functions (diagonal and off-diagonal ones) and the 
vertices, together with the resonance condition at each vertex, 


we end up with the following relations 

zqQ^zq^Q', (5.13) 

zpP^ZP'P', (5.14) 

ZkK^Zk'K', (5.15) 

zqQ = -ZkK - zpP, (5.16) 

Zglql + Zp\p\ = -ZKlzgq + Zpp\, (5.17) 

ZkZk' + ZqZq’ = 2zpZp' . (5.18) 


Eq. ( |5.18| l is a result of the two independent resonance con¬ 
ditions in Fig. 1^ and reduces the number of diagrams for the 
retarded self-energy significantly. As a result of Eq. ( |5.18| l, 
the product of all the corresponding z-factors must be iden¬ 
tical. For the diagonal self-energy, where zqZq = -1, this 
means that zpZp' - ZkZk' = -1 and consequently only loops 
consisting of two diagonal Green’s functions contribute to the 
diagonal self-energy. Consequently, the diagrammatic repre¬ 
sentation of the diagonal self-energy is still given by the loops 
shown in Fig. even in the presence of off-diagonal Green’s 
functions, i.e. 

“ ^ fp “ pISpSq-p + If + ^ISq+pSp 

+lp + llgfgppg) ■ (5-19) 


with the off-diagonal self-energies Fg and the Green s func- jjj order to obtain the off-diagonal self-energy F*, one has to 
tions gg = -zXQrgQ'g),g^ = -zXagag),/zg = -/(agaXg) and flip the sign of zq> -> 1, such that zqZq' = ZpZp> = ZkZk' = 
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1 in the corresponding loops. This means that the diagrams 
contain only off-diagonal terms. Tg is then obtained by the 
diagrams in Fig. but with all arrows from the right vertex 
flipped. We thus obtain 

r 

~ '^0 " ^ jp " P\hphQ-P + \P + q\hQ+php 


-\p + q\{hff hl^o). (5.20) 


The diagonal self-energy 2* diverges when the integral ( |5.19| l 
is performed with the bare Green’s functions. This hints that a 
non-trivial self-energy is generated on the diagonal to regulate 
the integral, which can be computed in self-consistent Born 
approximation as explained in previous sections. On the other 
hand, for off-diagonal Green’s functions - h'^ -0 (e.g. for 
the bare off-diagonal values), the off-diagonal self-energy is 
zero, as visible from Eq. ( |5.20| l. Consequently, off-diagonal 
self-energies are not generated in the retarded sector even in 
self-consistent Born approximation and we have 

r« = 0. (5.21) 

In the absence of off-diagonal self-energies in the retarded 
sector, we can directly apply the quasi-particle approxima¬ 
tion discussed in the previous sections and evaluate the self¬ 
energies and distribution function on-shell. Consequently, the 
intermediate kinetic equation is 

dTFg^^=o,T = q,oj=0,T + io-z^lo,=o,T^z (522) 

with the scalar, on-shell self-energy as discussed in Sec. 

mm 

For the on-shell Keldysh self-energy. 


^q,cj=0,T 





YK 

g,T 


(5.23) 


the diagrammatic rules from the previous section do not have 
to be modified. This immediately yields the di^onal, on- 
shell Keldysh self-energy according to Fig. For the 
off-diagonal, on-shell Keldysh self-energy T^ the correspond¬ 
ing diagrams are obtained by reversing the arrows associated 
to the vertices on the right in Fig. resulting in the diagram¬ 
matic representation of T^ shown in Fig. ^ The Keldysh 
Green’s functions are obtained via the parametrization used 
in Eq. (|5.9|), which yields for the off-diagonal elements 


h 


K 


len^io-l^nig^r 

(oj - u\q\)^ + 


(5.24) 


Evaluation of the off-diagonal diagrams and insertion into the 
kinetic equation completes the set of equations for a system 
including anomalous phonon densities. 

Finally, the off-diagonal retarded/advanced self-energies 
are exactly zero and the kinetic equation for the off-diagonal 
terms is 


dfing = 
-I- 


/ 

f 


2pq(q - p) 

o<p<q d-’^ + &l + 


Apq(q -I- p) 


+ cr‘' + cr‘ 


q+p 


{lUpWlg^p -mq{\ 
p+qf^q "I" ^^^p^^p+q 


+ Up + f^q-p^ 

- rngiip) .(5.25) 



FIG. 9. Diagrammatic representation of the off-diagonal Keldysh 
self-energy F^. Compared to the diagonal component, for the vertex 
of the left ingoing and outgoing lines have heen replaced. The di¬ 
agrams containing only retarded/advanced Green’s functions, which 
contributed to cr^, are absent because of the absence of off-diagonal 
retarded and advanced Green’s functions. 


Here, we have again used the transformed basis to elimi¬ 
nate the factor Vq in front of the integrals. The time evolu¬ 
tion for the off-diagonal terms depends on both the diagonal 
and off-diagonal densities, and the stationary solution of this 
equation is = 0 due to the uncompensated spontaneous 
decay term in the first integral. Together with Eqs. ( |3.15| ) 
and (EH', Eq. ( |5.25| l represents the complete set of equa¬ 
tions determining the time evolution of the phonon density 
and the self-energies of an interacting Luttinger Liquid for a 
system that has been initialized in an out-of-equilibrium state 

Hq + n^(u\q'\),mq + 0. 


VI. RELAXATION OF AN EXCITED THERMAL STATE 


In this section, we will analyze the relaxation dynamics of 
a nearly thermal initial state and compare it to the analytical 
results from the previous sections. We consider an initial state 
with the densities 

n^.T^o = «B(Ti = 2 m|^oI) + 5q 

- V a\e ’"2 , (6.1) 

- 1 

(g-gp)- 

mq^T=o = = 2aie ""2 . (6.2) 

A state of this form can be created by perturbing a thermal 
state with initial temperature T; = 0.5M|g'ol in coupling the op¬ 
erator cx dx(t> to a classical field with momentum qo, i.e. in 
a microscopic fermionic or bosonic model by a small density 
modulation with momentum g’o.^^For a specific simulation, 
we chose the parameters aj = 0.2, q '2 = A-q^ and express mo¬ 
mentum in units of qQ. The perturbation increases the energy 


0-5 - 
e ItoI 


£(t) = 



(6.3) 


of the system. As a result, the final state in the limit t ^ oa 
will be a thermal state with increased temperature Tf > T^. 
Since the kinetics is energy conserving, E{t) = E(t - 0) for 
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FIG. 10. Time evolved phonon density n,,. after initializing the sys¬ 
tem in a state described by Eqs. l |6.1[ >, l |6.2| l. The time t is expressed 
in units of . a) The deviation ^ of n„r from the final thermal 

density n^{Tf) divided by UefTf). decays to zero with a momen¬ 
tum dependent rate « -l.lcr^. For smaller momenta, this rate de¬ 
creases and finally becomes zero for ^ > 0. b) Anomalous density 

for different times t. The pinning of is the consequence of 
global particle number and current conservation of the present model. 


all T > 0, the temperature can be determined according to 
r u\q\ T^tP' 

from the system energy. In the present case, this leads to Tf - 
0.52M|q'o| and a final state of the system 


lim = naiTt)- (6.5) 

r-^oo 

The quantity of interest is the deviation of the time-dependent 
phonon density from the final phonon density in the limit 

T OO, 




( 6 . 6 ) 


In Fig. 10 we show j and m^ j for different time steps 
Tj, and we see in which way both quantities decay to zero 
momentum as a function of time and momentum. 

The evolution of this excited state to its new equilibrium 
proceeds in two different stages® For short times, the set of 
excited modes is much more strongly occupied than all other 



FIG. 11. Comparison of the decay rate jg = —dr log as a func¬ 
tion of momentum for different times t with the analytically esti¬ 
mated decay factor = 2.2(7^. The time is given in units of -K. 

The exact decay coincides very well with the analytical prediction. 
Only for 0, small numerical errors leading to drfig 0 lead to 
a divergence of the logarithm. 


modes, and the dominant effect is the decay of these modes 
into the continuum of non-excited ones. In this step, there 
is no back-action from the continuum, which acts as a bath 
for the excited states. As a consequence, the dynamics is de¬ 
scribed by the corresponding equilibrium linear response. Ac¬ 
cording to Eq. ( |4.11[ ), driq then follows the effective equation 
of motion 


log = -l.lo^. 


(6.7) 


As a result, for short times drig^r decays exponentially in time 
with a momentum dependent decay time - (l.lcr^) . 

; compare the numerical value of d-^ log [Srig ^ 


In Fig. 
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we ( 


with the analytical estimate 

log~-l-lo^,T ~ find a good 

agreement in the momentum region where 6nq deviates from 
zero. 

For larger times instead, many modes deviate only very 
weakly from the thermal occupation, but the system has still 
not found its equilibrium. In this case, back-action from the 
mode continuum can no longer be ignored. More precisely, 
due to energy and momentum conservation, this dynamics, 
which now is dominated by in- and out-scattering events on 
an equal footing, is very slow. It is no longer described by an 
equilibrium response theory, which is determined solely by 
the retarded self-energy. Instead the presence of dynamical 
slow modes is revealed, which are not captured by a simple 
exponential decay but have to be implemented as additional 
modes due to symmetrieJ^ and favor an algebraic decay in 
time, i.e. 

Srigir) ~ T~^ with some exponent (3. In Fig. 12 the two dif¬ 
ferent stages of the time evolution are visible in the numeri¬ 
cal simulation of the quasi-particle occupation. The algebraic 
decay for long times, and the exponential decay for shorter 
times, are clearly distinguished. In order to estimate the al- 
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FIG. 12. Decay of the deviation of the phonon density to 

its corresponding thermal value for fixed momenta {qi,q 2 ,qi) = 
(0.75,1,1.25)^0- a) Log-Log plot of the deviation Sn^ir) (blue dots: 
simulations; line: algebraic fit): For long times, the decay is de¬ 
scribed by a power law in time, i.e. 6n^{T) ~ with an exponent 
p X 0.58. b) Semi-logarithmic plot of the deviation SiiqiTy. For short 
times, the decay is described by an exponential, 6n^(T) ~ 
where is the corresponding thermal self-energy (blue dots: nu¬ 
merical simulation; red line: exponential fit, obtained from the self¬ 
energy). 


G = - 

, V = ^ 

r(4) = 



FIG. 13. Schematic diagrammatic representation of the self-energy E 
and the full three-point vertex V in terms of the full Green’s function 
G, bare three-point vertex V, and the full four-point vertex F*''*. 


and large momenta. 


gebraic exponent, we do not linearize the kinetic equation but 
instead take Eq. and insert on the right hand side the so¬ 
lution for the phonon occupations riqif) — iigiTf) ~ 

Due to the scaling of the thermal self-energy cr^ ~ qi, this 
yields the estimate p - |. In the numerical simulations pre¬ 
sented in Fig.[^ we find for the algebraic exponent p x 0.58. 

As mentioned already, this algebraic decay is a conse¬ 
quence of energy and momentum conservation in the dynam¬ 
ics, which leads to additional slow modes in the time evo¬ 
lution. The separation of the relaxation into two different 
time regimes, with first exponential decay according to bare 
phonon decay and then algebraic decay due to energy con¬ 
servation, has already been found in a recent work using an 
equilibrium formalism.!^ In this work, energy conservation is 
ensured by an effective, time-dependent temperature, which 
can, however, not account for local energy fluctuations. On 
the other hand, in the formalism presented in this manuscript, 
both spatial and temporal energy fluctuations are naturally in¬ 
corporated in the kinetic equation. The algebraic decay of the 
phonon occupations can be explained in terms of dynamical 
slow modes, resulting from conservation laws, i.e. symme¬ 
tries, in the system. The analytical exponent P - \ would 
belong to a system with exact momentum and energy con¬ 
servation with a classical distribution function n„ ~ jL xhe 
deviation of the numerical result p x 0.58 from the analytical 
one might be a result either of subleading corrections that will 
vanish on even larger time scaled^ or the fact that the system 
is not described by a classical distribution function ~ jL 
for the complete momentum range but by a Bose distribution, 
which seriously differs from the classical one at intermediate 


VII. DYSON-SCHWINGER EQUATIONS AND VERTEX 
CORRECTIONS 

In this section we apply Dyson-Schwinger equationd^iss] jq 
the interacting Luttinger Liquid and determine the self-energy 
and three-point vertex as a function of the phonon distribu¬ 
tion. We show that the vertex correction is always real and 
will be exactly zero for a zero temperature state. More gen¬ 
erally it leads to a modification vq ^ v'o (1 1)7 where J is a 

small (1 » I) dimensionless function with weak momentum 
dependence, whose precise form is determined by the time- 
dependent phonon density. Based on these findings, we con¬ 
jecture that also the corrections to the four-point and higher 
order vertices, which do not occur in the microscopic action, 
are small. As a result, the kinetic equation and the equation 
for the self-energies in self-consistent Born-approximation are 
modified according to the vertex correction, and we derive a 
coupled but closed set of equations. 

In general, Dyson-Schwinger equations (DSE) represent 
an exact hierarchy for one-particle irreducible (IPI) correla¬ 
tion fun ctions, generated by the effective action functional 
r[aQ,, .1321111] It relates the full vertices (e.g. the Green’s 
function, which is the inverse two-point vertex) to their bare, 
microscopic counterparts, which form the microscopic action 
S. The hierarchy built up by the DSE is in general infinite and 
therefore lacking an exact solution. The main goal is then to 
And a physically reasonable truncation of the effective action, 
for which the main physical effects are captured and in which 
the DSE can be solved. 

The effective action, as the generator of IPI correlation 
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functions, can 


be expandend according 


and for the full three-point vertex 


DO j 

r[a„,a/j] = ^ —^al...a,aai-aa„, (7.1) 

«=2 ”• 


Vafiy — + iS''^g^GsijVfij3,jGria-Va-ycG^v 




£ c (3) ^ 

2^ aSv'-'Sfi^ tiPyti^’l'" 


■(4) 


(7.5) 


where T^au-a„ the n-th order functional derivative of the ef¬ 
fective action. 


p(n) 


dUn-.-dai 


(7.2) 


i.e. the full n-point vertex function, and Oa are the fields with 
collective index a - (to, q, t, c/q). The inverse Green’s func¬ 
tion and the full three-point vertex are the second and third 
order functional derivatives, respectively. In terms of formu¬ 
las we have 


G^' = Go' - S = and V = (7.3) 

In the present case, the DSE relate the n-point vertex to the 
n + 1 point vertex of the theory. More specifically, we obtain 
for the full Green’s function 

= {^0% - iS%Gy,V,,.G.s, (7.4) 


Here, is the bare (microscopic) three-point vertex. A 
schematic diagrammatic representation of Eqs. ( |731 l, ( |74| 
and ( |7.5| l is depicted in Fig. 

We will show in the following that there are corrections to 
the bare three-point vertex, with the same scaling dimension 
as the bare three-point vertex itself. It is determined by a di¬ 
mensionless function I whose precise form is dictated by the 
time dependent phonon density This is in contrast to the 
two-point vertex, i.e. the inverse Green’s function, where the 
correction due to the cubic vertex is subleading but introduces 
an imaginary part and therefore is of great qualitative impor¬ 
tance. Here the correction is purely real, as the bare vertex 
itself, but on the other hand not subleading. It therefore can 
not be discarded directly without further discussion. 

For the ag^pUpOQ term, the vertex correction 6 V = V 
is illustrated in a diagrammatic representation in Fig. It is 
equivalent to all other vertices that incorporate a single quan¬ 
tum and two classical fields. The lowest order contribution, 
incorporating only bare vertices, reads (for q,p > 0) 


13 


q,p,p+q 


= -p yl\qpip + q)\ f ^ 

V8 Jk>0 


; k(q + k)(k - p) 


k+q 




k-p 


Hk- 


■rik 


^q+k ^k 


CTj + cr‘ 


k-p 


^k+^kyg 


k(p + k)(k - q) 
erf -I- erf 

k+p k-q 


Hk- 


Hk 


^p->rk ^k 


+ ^k 




(7.6) 


Here erf = -Im(2f^) > 0 are the on-shell self-energies ob¬ 
tained by the DSE in Fig. 13 The vertex correction due to 
cubic vertices is identical to zero for a zero temperature sys¬ 
tem (rik = 0 for all momenta). This is an exact statement for 
the interacting Luttinger Liquid and can be seen on the level of 
the diagrams in Fig. 14 For T - 0, G^ = sgn(fa;) (g* - G^^ 
and the individual diagrams cancel each other due to the pole 
structure of G* and G'^. In general, for a constant quasi¬ 
particle distribution function, the cubic vertex correction is 
exactly zero in arbitrary dimensions. 

In order to find a compact expression for the vertex cor¬ 
rection, we replace cr* ^ vqij^ and compare the integrand 
in E q. (|7.6| l with the expression for the self-energy in 
Eq. ( |3.18 1. We immediately see, that the vertex correction 
is linear in vq and the integral has scaling dimension zercPl 
In perturbation theory, this yields the vertex 


Here, Jo (|, n) is a dimensionless function of the ratio pjq and 
the phonon density n, which is determined by the integral in 
Eq. (1^. The scaling behavior of the one-loop vertex correc¬ 
tion suggests the parametrization of the full vertex according 
to 


yq,P,q+P = VQ y[\pq(pVq)\ 



(7.8) 


y(lst) 

''q,p,q+p 


Vo 


^Jlpqip + q')| 1 + Jo 



c(3) 

*4 g,p,p+q 


1+Jo 



where the functional J encodes the full vertex correction. Ac- 
(7.7) cording to the DSE in Fig.[T^ it is determined via 
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I{x, n) 
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0 -? H- 0 -? erf -1- o"? 
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(7.9) 
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FIG. 15. Vertex correction I(x,n^(T)) for an infinite temperature 
state. In the limit T —» oo, the temperature drops out in Eq. .9) , and 
the vertex correction becomes temperature independent. Due to the 
invariance of I under x —> f the plot is restricted to 0 < x < 1. The 
dependence of J on x is weak, especially for x ~ 1, where it takes its 

maximum J(x, ng) ^ Ima ~ 0.093. 


FIG. 14. Diagrams contributing to the vertex correction of the three- 
point vertex illustrated for the particular example of an outgo¬ 

ing quantum field and two incoming classical fields, i.e. oc 
In total there are twelve distinct diagrams contributing to the vertex 
correction, six are depicted above and six further can be found by 
interchanging the ingoing momenta P 2. 


where ^ ^ and the integral is dimensionless. In Eq. 

we have already exploited the fact that the ingoing momenta 
of a vertex can be exchanged without modifying the vertex 
itself, and consequently the integral is invariant under q ^ p. 
This is equivalent to J(x, n) = n). The self-energy in the 

DSE approach is 



(7.10) 



qp{q - p) 


cr" -H cr; 


q-p 



In the same way the kinetic equation can be modified to in¬ 
corporate the vertex correction as well, and we find a set of 
coupled equations 


' drll 
(T* 

I 


= Tin, 0^,1). 


(7.11) 


They can be solved numerically according to the procedure 
described in Eig.|^ Including the vertex correction, the second 
step of the iteration additionally includes the self-consistent 
determination of I. 

We will now give an estimate of the order of the vertex cor¬ 
rection for the case for which it is non-zero to estimate its 


impact on the dynamics and the kinetic equation. In the limit 
T —» oo the relevant phonon density is riq ^ and the self¬ 


Consequently, 


energies have the thermal form of oc 

the factor ^ drops out and I (x, n) does no longer depend on 
temperature. In this case, 1(0,n) 0.012, J(l,n) ^ 0.09 and 
-^(0, n) < I(x,n) < I(\,n). As can be seen in Eig. 15 the cor¬ 
rection is small, with a weak momentum dependence. Conse¬ 
quently, the self-energy is only negligiably modified if instead 
of the full three-point vertex in Eq. the bare value V 

is used. This is precisely the self-consistent Born approxima¬ 
tion that we used to determine the self-energies and the kinetic 
equation for the interacting Luttinger Liquid. 


VIII. CONCLUSION 

In this article, we used non-equilibrium field theory, in par¬ 
ticular kinetic and Dyson-Schwinger equations, to determine 
the kinetics and non-equilibrium dynamics of resonantly in¬ 
teracting Luttinger Liquids. Exploiting the fact that the inter¬ 
actions lead to dressed but still well defined phonons, which 
enables a separation of timescales into slow forward and fast 
relative dynamics, we applied the Wigner and quasi-particle 
approximation and derived a closed set of simple yet power¬ 
ful equations for the normal and anomalous phonon density, 
the phonon self-energy and vertex correction. These equa¬ 
tions determine the dynamics of an interacting Luttinger Liq¬ 
uid initialized in a Gaussian (non-) equilibrium state. The re¬ 
sulting equations show strong aspects of universality, on the 
one hand being independent of any UV-scale, in particular in¬ 
dependent of the Luttinger cutoff. On the other hand, after a 
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proper rescaling of the forward time, all microscopic parame¬ 
ters entering the Hamiltonian can be eliminated and the only 
microscopic information entering the dynamical equations is 
the initial phonon density. We further used our approach to 
analytically determine the relaxation rate of a thermally ex¬ 
cited state. For this dynamics, we found an initial exponential 
decay correspo nding to pre vious results computed from linear 
response theory.EEHSIlSIZl] However, for longer times, the de¬ 
cay of the excitations follows a power law in time, revealing 
the presence of dynamical slow modes due to energy conser¬ 
vation. These latter modes are not contained in a plain equilib¬ 
rium linear response theory but have to be build in by ha nd on 
the basis of symmetry arguments and conservation laws.l^^HS] 
Here, the dynamics based on the kinetic equation approach re¬ 
veals the presence of these modes due to the algebraic decay at 
long times without any further modification, which shows the 
strength of our approach in a simple yet nontrivial example. 

The results of this work can be used in order to de¬ 
termine the kinetics and non-equilibrium dynamics of one¬ 
dimensional interacting quantum fluids prepared in a non- 


thermal initial state, which might occur as a consequence of 
a quantum quench or a sudden external perturbation. On the 
other hand, it paves the way to compute the dynamics of quan¬ 
tum fluids subject to drive and dissipation, and to determine 
the dynamics towards the steady state of an excited closed 
system. Of special interest and a strength of our approach is 
the treatment of long time dynamics in non-equilibrium sys¬ 
tems, which are neither reachable by present numerical pro¬ 
cedures nor by analytical approaches based on perturbation 
theory. Both examples belong to the uprising held of one¬ 
dimensional quantum fluids out-of-equilibrium and we leave 
their discussion open for future work. 
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